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We study ultranarrow 5S]/ 2 - 4/ 5 / 2 transitions in Nd 13+ and Sm 15+ and demonstrate that they 
lie in the optical region. The transitions are insensitive to external perturbations. At the same time 
they are sensitive to the variation of the fine structure constant a. The fractional accuracy of the 
frequency of the transitions can be smaller than 10 -19 , which may provide a basis for atomic clocks 
of superb accuracy. Sensitivity to the variation of a approaches 10 -20 per year. 
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Building accurate atomic clocks is important for both 
technical applications and fundamental physics. Cesium 
primary frequency standard which is currently used to 
define the SI units of time and length has fractional ac- 
curacy of the order 10 -16 [l|. Frequency standards based 
on optical transitions in neutral atoms trapped in opti- 
cal lattice aim at fractional accuracy of 10- 18 @. Fur- 
ther progress is possible with clocks using nuclear optical 
transition Q , or clocks using optical transitions in highly 
charged ions 

One of the important applications of precise atomic 
clocks in fundamental physics is the search for time vari- 
ation of the fine structure constant. The possibility for 
fundamental constants to change in space or time is sug- 
gested by theories unifying gravity with other interac- 
tions (see, e.g. 0,@])- The analysis of quasar absorption 
spectra shows that there might be a direction in space 
along which the fine structure constant a (a = e 2 /he) 
is changing at constant rate over cosmological distances 
(a-dipole) @. Since the Earth moves in the framework 
of the a-dipole, the space variation of a should mani- 
fest itself in terrestrial studies as changing of a in time 
at approximate rate of a/a = lO^^yr -1 10]. The 
best current limit on time-variation of a obtained from 
comparing Al + and Hg + clocks is significantly larger, 
a/a = (-1.6 ± 2.3) x lO^V -1 [ll|. At least an or- 
der of magnitude improvement is needed to verify the 
a-dipole hypothesis in terrestrial studies. Two or more 
clocks with different dependence of clock frequencies on 
the fine structure constant and operating on the same 
level of accuracy (10~ 18 or better) should be monitored 
over significant time interval (from few months to few 
years) to reveal any variation of a in time. If a is chang- 
ing, than the readings of the clocks would shift in time 
against each other. 

It was suggested in Ref. [|[ to use highly charged ions 
(HCI) for this purpose. HCIs are less sensitive to exter- 
nal perturbations than neutral atoms (or singly-charged 
ions) due to their more compact size. Notice that many of 
HCIs have optical transitions where standard laser tech- 
niques can be used and where accuracy of the frequency 



measurements can be possibly even better than in optical 
standards using neutral or single-ionized atoms. A num- 
ber of atomic transitions in HCI were studied in earlier 
works [U, [H, [l2T - [l4j focusing mostly on their sensitivity to 
the variation of the fine structure constant. A detailed 
study of the effect of external perturbations on clock fre- 
quencies of some promising HCI systems was recently 
reported in Ref. @ . There it was shown that the relative 
accuracy on the level of 10 -19 might be possible for HCIs 
with the 4/ 12 configurations of valence electrons. These 
systems are not sensitive to the variation of a. However, 
they can be used as a reference point against which the 
frequencies of the transitions sensitive to the variation of 
a are monitored. 

In present work we study other promising systems 
which arc immune to external perturbations while be- 
ing sensitive to the variation of the fine structure con- 
stant. This makes them excellent candidates for terres- 
trial studies of the time variation of a. We focus on the 
5s - 4/ transitions in Nd 13+ and Sm 15+ ions. Frequencies 
of these transitions are in the optical region due to the 5s 
- 4/ level crossing which happens for the isoelectronic se- 
quence of silver atom near Z = 61 (promethium ion) [Til ]. 
The first excited states in these ions (4/ 5 / 2 in Nd 13+ and 
the 5s in Sm 15+ ) are metastable states because they can 
only decay to the ground state via the E3 transition. 
The states are not sensitive to external perturbations. 
On the other hand, the s — f transitions are the transi- 
tions which are most sensitive to the variation of the fine 
structure constant [ll| . Due to the inversion of the states 
from Nd 13+ to Sm 15+ , corresponding frequencies would 
move in time in opposite direction if a changes. This 
further enhances the sensitivity. If the frequencies of the 
5s — 4/5/2 transition in Nd 13+ and 4/ 5 / 2 — 5s transition 
in Sm 15+ are monitored simultaneously for extended pe- 
riod of time the sensitivity of the measurements can be 
on the level 8a /a ~ 10~ 19 per year or better. 

Note that similar transition in Pm 14+ is even more 
sensitive to variation of a. However, in this work we 
consider only isotopes with zero nuclear spin to suppress 
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the second-order Zeeman shift. The promethium atom 
has no stable isotopes and no long-living isotopes with 
zero nuclear spin. Therefore we limit ourself in this work 
to the Nd 13+ and Sm 15+ ions. 

The sensitivity of the atomic transitions to the vari- 
ation of the fine structure constant can be revealed by 
varying the value of a in computer codes. 

We use the correlation potential method [16| to per- 
form the calculations. The Nd 13+ and Sm 15+ ions have 
a single valence electron above the Pd-like closed-shell 
core. Therefore, it is convenient to use the V N_1 ap- 
proximation in which relativistic Hartree-Fock calcula- 
tions are first done for the closed-shell core and states of 
valence electron are calculated in the self-consistent field 
of the frozen core. Correlations between valence and core 
electrons are included by constructing the second-order 
correlations potential E and solving the Hartree-Fock- 
like equations with an extra operator £ for the states of 
external electron 



(Hq + E — e v )il> v = 0. 



(1) 



Here Hq is the relativistic Hartree-Fock Hamiltonian and 
index v labels valence orbitals. Solving ((TJ) gives energies 
and wave functions for different states of the valence elec- 
tron which include correlations (the so called Brueckner 
orbitals). Many-body perturbation theory and B-spline 
basis states [13] are used to calculate X. 

The dependence of the atomic frequencies on the fine 
structure constant a the frequencies may be presented as 



uj(x) — ojq + qx, 



(2) 



where x = (a/ao) 2 — 1, a is the current value of the 
fine structure constant, and ao is the value of the fine 
structure constant at some fixed moment of time, say the 
beginning of the observations. Note that we use atomic 
units which means that the unit of energy is fixed and 
does not vary with a. Therefore, equation ([2} gives com- 
prehensive description of the frequency dependence on a. 
In the end, only change of the dimensionless values (e.g., 
ratio of two frequencies) can be studied. Therefore, the 
actual choice of units cannot affect the results. 

The sensitivity coefficient q in © is calculated by run- 
ning the computer code with different values of a and 
taking numerical derivative, 



w(0.01) -w(-0.01) 
O02 ' 



(3) 



It follows from © that the relative change in frequency 
is related to the relative change in a by 



8lo 2q 5a 



id a 



(4) 



where K = 2q/u is the enhancement factor. 

We chose the Nd 13+ and Sm 13+ ions because the 5s - 
4/ transition for these ions lie in the optical (UV) region. 
Table Q] shows calculated energy levels of the 4/ and 5s 



TABLE I: The 5s and 4/ energy levels and sensitivity coef- 
ficients q for Nd 13+ and Sm 15+ (cm -1 ). K = 2q/ui is the 



enhancement factor. 


Ion 


State 


Energy 


1 


K 


Nd i3 + 


5Si/ 2 













4/5/2 


58897 


106000 


3.4 




4/7/2 


63613 


110200 


3.5 


Sm 15 + 


4/5/2 













4/7/2 


6806 


6300 


1.9 




5Si/ 2 


55675 


-136000 


-4.9 



states of the ions together with the sensitivity coefficients 
q (see ((2]) and (fTTj) ) and the enhancement factor K. All 
other states of these ions lie very high in the spectrum, 
far outside of the optical region. The 5s - 4/ 5 / 2 transition 
is very narrow for both ions. The lowest-order transitions 
between the 5s and 4/ 5 / 2 states are the M2 and E3 tran- 
sitions. The M2 transition is very weak since it vanishes 
in the non-relativistic limit. The width of the lines is de- 
termined by the 4/ 5 / 2 - 5s E3 transition for Nd 13+ and 
by the 5s - 4/ 5 / 2 and 5s - 4/ 7 / 2 E3 transitions for Sm 15+ 
ion. The E3 decay rate is given by (we use atomic units: 
H = 1, m P = 1, lei = 1) 



T h = 0.001693aV, 



{a\\E3\\b) 2 
2J b + l 



(5) 



Here a is the ground state and b is the excited metastablc 
state, a = 1/137.36 is the fine structure constant, Lu a b 
is the frequency of the transition. Estimates show that 



for Nd 13 + T 



4/5/2 



1.1 x lCT^ 3 a.u. = 7.4 x lCT a Hz, 



l/Q = T ih/2 /oj w 4.8 x 1CT 23 (Q is the quolity factor), 
and for Sm 15 + T 5s w 4.3 x lfT 23 a.u. = 2.8 x 1(T 7 Hz, 
l/Q = r 5s /u w 1.7 x 10~ 22 (see also TablelH]). 

Note, that the 4/ - 5s level crossing [1J] happens at 
about Z = 61 (promethium ion) leading to inversion of 
the levels order in samarium ion. As a results, the q- 
coefficients and enhancement factors have different sign 



for Nd 13+ and Sm 1 



It may be beneficial to mea- 



sure the frequency of the 4/ 5 / 2 — 5s transition in Sm 15+ 
(A = 170 nm) against the frequency of the 5s— 4/ 5 / 2 tran- 
sition in Nd 13+ (A = 180 nm). This would bring extra 
enhancement to the sensitivity of the clock frequencies to 
the variation of the fine structure constant. Using data 
from Table Q] one can get 



qi-q 2 = 242000 cm" 
6(uii/l} 2 ) ^ S5 f^_ 
(wi/w 2 ) 



(6) 



Here uii is the clock frequency for Nd 13+ and w 2 is the 
clock frequency for Sm 15+ . Eq. (J5]) shows that if a drifts, 
the ratio of the clock frequencies changes about eight 
times faster. This means that if the ratio is monitored 
for a year with the accuracy of 10 -19 , the sensitivity of 
the measurements to the variation of a would be close to 
10 -20 per year. 
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TABLE II: Parameters of the clock transitions for Nd 13+ and 
Sm 15+ . Wavelength of the transition (A), radiative width 
of the excited clock state (r), static dipole polarizabilities 
(a(0)), magnetic dipole polarizabilities (7c), and magnetic 
dipole hyperfine structure constant (A, gi = fJ>/I)- Numbers 
in square brackets represent powers of 10. 





Nd 13 


f 


C 15 

bm 


4- 




Ground 


Excited 


Ground 


Excited 




state 


state 


state 


state 




5s 


4/5/2 


4/5/2 


5s 


A [nm] 


170 




180 




r [Hz] 





7.4 [-8] 





2.8[-7] 


a(0) [ag] 


1.1048 


0.3701 


0.2766 


0.8366 


7c [a.u.] 





0.144 


1.344 





Ajgi [MHz] 


123000 


1070 


1400 


154000 



Below we consider various systematic effects which af- 
fect the clock transition frequency. As discussed in [f| the 
clock HCI is assumed to be trapped and sympathetically 
cooled. 

a. Black body radiation shift - The frequencies of 
the clock transitions might be affected by the black body 
radiation (BBR) shift. The BBR frequency shift at an 
ambient temperature T can be expressed as 



8lj 




\ 4 Aa 






(f 




( 8 31. 9 X) 




\ 1 o / 


/ 2oj 





(7) 



where To = 300 K, Aa is the difference in the values of 
the static dipole polarizabilities of the clock states. The 
polarizability a v (0) for the state v is given by 



a v (0) 



3(2j„ 



7>£ 



(w||D||n) : 



(8) 



where D = J2i er i i s the electric dipole operator and 
summation goes over complete set of states. The re- 
sults of calculations, which include Brueckner-type corre- 
lations and core polarization effects, are presented in Ta- 
ble |TTJ Using (J7]) and the data from the table one can get 
at room temperature Slu/lo = —3.6 x 10~ 18 for Nd 13+ and 
Suj/uj = 2.9 x 10~ 18 for Sm 15+ . For cryogenic Paul trap 
operating at the temperature of liquid helium (~ AK) 
the fractional BBR shift is 5ui/lu ~ 10~ 25 for both ions. 

b. Zeeman shift - Clock frequencies are affected by 
magnetic fields. The first-order Zeeman shift can be elim- 
inated by averaging the measurements over two virtual 
clock transitions with opposite g-factors. Uncontrollable 
second-oder AC Zeeman shift arises due to misbalances 
of currents in ion traps. It can be evaluated as 



1 



2j v + 1 ^ Ei - E c 



B 



(9) 



The second-order Zeeman shift is negligibly small for the 
5s states since there are only strongly forbidden Ml tran- 
sitions in ^ . The shift for the 4/ 5 / 2 states is dominated 



by the 4/7/2 - 4/ 5 / 2 transition within the fine structure 
doublet. The calculated values of 7 C for these states are 
presented in Table [ill Using the value of AC magnetic 
field B = 5 x 10" 8 T measured in Al+/Be+ trap [TJ| 
leads to Suj/lj ~ 10~ 26 . 

c. Electric quadrupole shift- Clock frequencies can 
be affected by coupling of atomic quadrupole moments 
to the gradients of trapping electric field. Obviously, this 
does not affect the 5s state due it's vanishing quadrupole 
moment, but can shift the 4/ 5 / 2 state. It was suggested 
in Q to use hyperfine structure (Ms) to suppress the 
quadrupole shift. However, using isotopes with non- 
zero nuclear spin leads to enhanced second-order Zeeman 
shift. This is because small hfs intervals enter energy de- 
nominators in ^ increasing the shift by several orders of 
magnitude compared to isotopes with zero nuclear spin. 
For this reason in this work we propose a different method 
to suppress electric quadrupole shift. It uses a linear 
combination of transition frequencies between states with 
different projection of total angular momentum J. 

The quadrupole shift for a state with total angular 
momentum J and its projection J z — M is given by 



SE. 



JM 



3M 2 - J( J + 1) dE z 



2J(2J- 1) 



0-. 



cjmQ 



dE z 
dz 



(10) 



where Q is the electric quadrupole moment. The depen- 
dence of the shift on the projection M is in the prefactor 
only. The quadrupole moment Q of the atomic state 
is defined as twice the expectation value of the electric 
quadrupole operator (E2) in the stretched state 



Q = 2(nJM = J\E2\nJM = J). 



(11) 



Taking two transitions to states with different val- 
ues of M and writing the frequencies as wm = u a + 
cjMQ(dE z /dz) one can get 



UJ 



UJM — WM'CJm/cJM' 
1 — Cjm/Cjw 



(12) 



Here loq is the frequency of the transition at zero electric 
and magnetic fields. The expression (fl"2j) does not depend 
on the quadrupole moment nor the gradient of electric 
field. For better accuracy it is important to have cjm 
and cjm> as different as possible. Using €5/2,1/2 = —0.4 
and C5/2.5/2 = 0.5 leads to 



UJ 



h/2 



0.8a; 



5/2 



(13) 



d. Other perturbations - The performance of ion 
clocks can be affected by many other systematic ef- 
fects, such as density of the background gases, Doppler 
(motion-induced) effects, gravity, etc. Consideration sim- 
ilar to what was done in our previous workQ show that 
none of the corresponding fractional frequency shift is 
below the value of 10 -19 . 
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To summarize, the Nd 13+ and Sm 15+ HCIs may offer 
an intriguing possibility for developing clockwork of un- 
precedented accuracy that is highly sensitive to variation 
of the fine-structure constant. The 10 -19 fractional ac- 



curacy matches projected accuracy of the 229 Th nuclear 
clock 3]. but without complications of radioactivity. 
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